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INTRODUCTION 

MIMIC is a programming system written for a digital 
computer which, from the standpoint of the user, seems to 
make the sequential machine function like the parallel 
analog computer and at the same time virtually eliminates 
the problems of time and amplitude scaling. It does this 
by providing a language which is based on the idea of 
interconnecting basic computing elements such as adders, 
integrators, switches, etc., as in the wiring of a 
patch board for a real analog computer. In this respect 
MIMIC is like MIDAS and PACTOLUS but has a larger set of 
computing elements, is faster, more efficient and easier 
to use for it allows FORTRAN-like expressions, MIMIC 
language sub-programs and logical control of program 
execution. 

This manual presupposes a limited acquaintance 
with analog techniques on the part of the reader. 



2. 



II. DEFINITIONS 

The following definitions will be used throughout 
this manual. 

1. Variable Name - A group of 1 - 6 letters or digits of 
which one must be a letter. 

e.g. AB, IB2 , 24C2 

The following symbols have special meanings to the 
MIMIC processor: 

T, DT, DTMAX, DTMIN, TRUE, FALSE. 

2" Literals - There are two types of literals: 

(a) Numeric constants which must conform to the rules 

(i) a decimal point 
(ii) at least one digit 

/j4^\ -.4- m^^-i- ,-.-:-,,. ^.'u-.-u.-t^-i-i-^-^-^ 

\j.ij.; a. l. niuo l OJ.A ^naia^ucio 

e.g. 47., 63.0, 4.5278. 

(b) Logical constants TRUE, FALSE. 

3. Function - A three letter mnemonic code for a MIMIC 
computing element. All MIMIC computing elements are 
specified in Table 1, Section III. 

e.g. LOG, ADD. 

4. Expression - A MIMIC expression may be: 

(i) a literal 

or (ii) a variable name 

or (iii) a MIMIC function with its arguments enclosed 
in parentheses and separated by commas. 

or (iv) any algebraically meaningful combination of 
these formed using the arithmetic operators 
+ - */** and parentheses. 

NOTE : Multiplication cannot be implied, 

e.g. A* (B + C) 

is legal but A (B + C) is not. 
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e.g. 28.0 
-IBX 
LOG(X,-2.) 

1BX+28 . 0*SIN (T* (THETA1-THETA2 ) /360 . ) 
NOTE : The arguments of MIMIC functions may be expressions. 

5. Statement - A MIMIC statement consists of two things: 
an expression and a variable name. The effect of the 
statement is to assign the result of evaluating the 
expression to the variable named. 

e.g. X ADD(Y,Z,-4. ,T*2 . 0+NEG (IBX) ) 

The result of carrying out the operations specified 
by the expression will become the value of the 
result variable named X. 

Exceptions: CON, PAR, RSP, OUT, HDR, FIN, END state- 
ments do not require a result variable. 

6. Program - A MIMIC program consists of a set of MIMIC 
statements, of which at least one must be a FIN 
statement, followed by an END statement. 

7. Job - A MIMIC job consists of a MIMIC program punched 
on cards and possibly followed by data cards. 

**• Processor - The MIMIC processor is a set of routines 
written in FORTRAN- IV and COMPASS for the 3300 which 
assembles and executes MIMIC programs. 

9. Special Variables - The MIMIC processor reserves the 
following symbols for the special purposes listed: 

T - the independent variable (considered possibly 
as time) 

DT - the amount T changes between printouts 
of values 

DTMAX - maximum integration stepsize allowed 

DTMIN - minimum integration stepsize allowed 

10. Program Run - The generation of the solution (to a 
problem programmed in MIMIC) from T = O until 
termination by a FIN statement will be called a run. 

The processor performs successive runs until it exhausts 
the supply of data cards provided for all PAR or PFN 
statements in the program. 



III. TABLE OF MIMIC COMPUTING ELEMENTS 



TABLE 1. LIST OF FUNCTIONS 



FUNCTION 


CODE 


1 INPUT (S) 




VALUE/REMARKS j 


1. 


ARITHMETIC FUNCTIONS 












Addition ADD 


A,B(,C,D 


,E 


rF)* 


R = A+B(+C+D+E+F) 




SUM 












Subtraction SUB 


A,B 






R = A-B 




Multiplication MPY 


A,B(,C,D 


,E 


rF) 


R = A*B(*C*D*E*F) 




Divide DIV 


A,B 






R = A/B 




Multiply and Add MAD 


A,B,C(,D 


r E 


rF) 


R = A*B+C(*D+E*F) 




Negation NEC 


A 






R = -A 




i »jw> •_» v-r _i_ vi «-. v_ vuxuc riJDCJ 
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Equality EQL 


A 






R = A 


2. 


ELEMENTARY TRANSCENDENTALS 












Square Root SQR 


A(>0) 






R =JK 




Sine SIN 


A(rads. ) 






R = sin A 




Cosine COS 


A(rads. ) 






R = cos A 




Arctangent ATN 


A(,B) 






R = tan" 1 (A/B) . 
If B is not 
specified it 
is assumed to 
be +1. 




Exponential EXP 


A(,B) 






R = B A . If B is 
not specified 
B = e is 
assumed. 




Logarithm LOG 


A(,B) 






R = log B A. If B 
is not spec- 






















ified B = e 












is assumed. 


* 


Operands enclosed in parentheses 


need not be 


spec 


if ied. 


1. 


B,C,D may be any lo 


gical-valued « 


expression 
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FUNCTION 



CODE 



INPUT (S) 



VALUE/REMARKS 



LOGICAL FUNCTIONS 
Function Switch 



Logical Switch 



And 



Exclusive Or 



FSW 


A,B,C.D 


R 


= B A < 
= C A = 




— 




= D A > 


LSW 


A,B,C X 


R 


= B if A TRUE 
= C if A FALSE 


AND 


A,B(,C,D,E,F)^ 


R 


=/TRUE if A and 
IB (and C and D 
land E and F) 
jhave value 
[TRUE* 

\FALSE, other- 
wise 


EOR 


A^B 1 


R 


=/TRUE if A and 
IB are differ- 
ent. 
/FALSE if A and 



Inclusive Or 



IOR 



Complement 



COM 
NOT 



A,B(,C,D,E,F) 



A - 



R = TRUE if A or 
B (or C or D 
or E or F) ha^ve 
value TRUE. 
Otherwise R = 
FALSE . 

R = (FALSE if A = 
I TRUE. 

[TRUE if A = 
i FALSE. 



INPUT/OUTPUT FUNCTIONS 

Name Constants CON 

Name Parameters PAR 



A(,B,C,D,E,F) Enters constant 

names 
2 
A(,B,C,D,E,F) Enters parameter 

names 



1. A,B,C,D,E,F must be logical-valued expressions 

2. Arguments must be variable names 

3. Argument must be a numeric literal 

4. SW is a control constant. If SW = 0., the data points will be 
plotted. If SW = 1., straight lines will be drawn between 
points. If SW=2., a quadratic will be fitted to the points. 
SCL is the multiplier to expand the scale in the X direction. 
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FUNCTION 


CODE 


1 INPUT (S) 


VALUE/REMARKS 


Name Function 


CFN 


A 3 


The name of the 


(Constant) 






array is entered 
in the result 
column. A = 
number of pairs 
or triples of 
points . 


Name Function 


PFN 


A 3 


The name of the 


(Parameter) 






array is entered 



Print OutDut 



Print Headers 



Plot 



OUT 
PRI 

HDR 

TTTTI TV 

PLO 



(A,B,C,D,E,F) 
(A,B,C,D,E,F) 



in the result 
column. A = 
number of pairs 
or triples of 
points. 

Print -A,B,C,D,E, 
F every DT units 
of T. 

Print heading 
names given in 
A,B,C,D,E,F 



SW,SCL,A(,B,C,D) Supply A,B,C,D to 

subroutine PLO 
every DT units 
of T. 



SUBPROGRAMS 

Begin Subprogram 
End Subprogram 



BSP 
ESP 



A(,B,C,D,E,F) 
A(,B,C,D,E,F) 



Call Subprogram 
Return Subprogram 



CSP 
RSP 



A(,B,C,D,E,F) 
A(,B,C,D,E,F) 



The subprogram 
name appears in 
the result column 
of the BSP and 
ESP cards. The 
inputs are named 
on the BSP card 
and the outputs 
on the ESP card. 

The name of the 
called subprogram 
is given in the 
result field of 
the CSP card. The 
RSP card must 
immediately follow 
the CSP card. 
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FUNCTION 



CODE 



INPUT (S) 



VALUE/REMARKS 



6. SPECIAL FUNCTIONS 



Integration 



INT 



A f B(,C,D) 



R = B + / AdT 

C,D must be 
logical val- 
ued expres- 
sions. 



\ c 

D \ 


TRUE 


FALSE 


TRUE 


OPERATE 


HOLD 


FALSE 


RESET 


OPERATE 

.,, ,i 



Limit Integrators 



First Order 
Transfer Function 



Limiter 



LIN 



FTR 



A,B,C,D 



A,B 



Dead Space 
Time Delay 



LIM 



DSP 



?DL 



A,B,C 



A,B,C 



A,B(,C) 



R = B < C 
= A C <, B <. D 
= B > D 

R = L™ 1 [A(s)/ 
(Bs+1) ] where A 
is the variable 
operated on by 
l/(Bs+l) 

R = B A < B 
= A B ^ A <. C 
= C A > C 



R = A - B 


A 


< B 


= B 


^ A 


<- C 


= A - C 


A 


> c 


R = A(T-B) 


. C 


is 


the number 


Of 




points of A to 


be 


stored and 


must be 


a literal < 


Dr con- 


stant. 






If C 


is not 


specified, 


C = 


100 


is assumed 


. B 


may 


be variable 


2. If 


T <L B, R = 


A(0) 


. 



& . 



FUNCTION 



CODE 



INPUT (S) 



VALUE/REMARKS 



Function 



Implicit Function 



Maximum 



Minimum 



Random No. Generator" 



FUN A,B 

A,B,C 

IMP A , B 



MAX A,B(,C,D,E,F) 
MIN A,B(,C,D,E,F) 
RNG A,B,C 6 



Random No. Generator" 



RNU 



Derivative 



_ _ _6 

A,B,C~ 



DER A f B,C 



R = A(B) 
R = A(B,C) 

R = A 

where |A~B (A) I 
£5x10 6 |A| 

R = max(A,B(,C,D, 
E,F)) 

R = min(A,B(,C,D, 
E f F)) 

R = random sample 
from a gaussian 
distribution 
with mean = A 
and standard de- 
viation B. C is 
a starting num- 
ber. 

R = random sample 
from a uniform 
distribution 
with lower limit 
A and upper limit 
B. C is a start- 
ing number . 



R = dA/dB. 
at T = 0. 



R = C 



HYBRID FUNCTIONS 

Monostable 
Multivibrator 



Track and Store 



MMV A,B 



TAS 



A,B,C 



R set TRUE when A 
TRUE and stays 
TRUE for B units 
of T after A 
goes FALSE. 
R = A at T = 0. 

R = A when B TRUE 
= R when B 
FALSE. 

(R = previous R) 
P 
R = C at T = 0. 
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FUNCTION 



CODE 



INPUT (S) 



Flip-Flop 



FLF 



A , B , C 



Zero Order Hold 



ZOH A,B 



8. CONTROL FUNCTIONS 



Go to Next Case 



End of Program 



FIN A,B 
END 



VALUE/REMARKS 



R =/ TRUE if A = 
TRUE . 

TRUE if B = 
FALSE and R = 
TRUE. P 
FALSE other- 
wise. 

R = C at T = 0. 

A is sampled every 
B units of T. R - 
held sampled value. 
R = A at T =0. 



Go to next run 

yvii^ii in. ^— j_» . 

Signifies end of 
MIMIC program and 
beginning of MIMIC 
data. 



5. 



6. 
7. 
8. 



Different starting numbers will produce independent sequences, 

hence many random number generators may be used in the same 

program. 

C should be a literal, constant or parameter. 

A,B,C must be logical valued expressions. 

R = A V(B'AR ) . 
P 
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IV. EXAMPLE 

The method of forming expressions and statements defined 
in section II means that- unlike MIDAS and PACTOLUS , only those 
variables of particular interest need be named in the result 
portion of a MIMIC statement. Thus a MIMIC program could be 
written down from the analytic form of the system being solved 
without the necessity of a block diagram. 

Suppose we wish to solve the problem 

x + x + x = with x(0) = 2, x(0) = 

A block (circuit) diagram representation of this problem 
is obtained by writing the highest derivative in terms of 
lower order derivatives, 



eg 
Thus we have 



x = 



-x -x 



ADD 



-x 



* NEG 



x 



jL 



3 INT- 



x 



I 



INT 



-& 



A MIMIC program which explicitly reflects the inter- 
connection of each basic computing element is: 



NEGDX2 


ADD(X,DX1) 


DX2 


NEG(NEGDX2) 


DX1 


INT(DX2,2.) 


X 


INT(DX1, 0.) 




0UT(X,DX1) 




FIN(T f 2.5) 




END 



The OUT and FIN statements cause the processor to print 
out a table of values of x, x for values of T = 0, .1, .2,.. 
2.5. 
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A simpler way of writing this program is 

DX1 INT(-(X+DX1) ,2.) 
X INT(DX1,0.) 

0UT(X,DX1) 

FIN(T,2.5) 

END 

V. ARBITRARY FUNCTIONS 

MIMIC allows arbitrary functions of 1 or 2 independent 
variables to be specified by tabulated data, i.e. by an array 
of values of independent variable (s) and corresponding values 
of the dependent variable. 

These functions are declared as such by including a CFN 

-1-i-nnN.T ^j.-.4-^™^-«-»- ■; « 4-Vn-> r>vArrram -F/~»v o^r«h funnt-ion rrpmPTa tor 

desired and are used in calculations by means of a FUN state- 
ment . 

For example the statement 

F CFN (5.) 

has the effect of defining a function F to be specified by 
5 sets of points which will be punched on data cards. Then 
the function F could be used, for example in the following 
statement 

X FUN(F,Z) 
i.e. the value of X will become F(Z) . 
Similarly 

X FUN(F,U,V) 
if F were a function of 2 independent variables. 

The processor uses linear interpolation for values of 
the independent variable (s) between tabulated points. The 
function is assumed continuous and to have derivatives equal 
to zero for all points outside the tabulated range. 
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The number of data points as well as the number of 
arbitrary functions is limited only by the computer memory 
available for programs. 

VI. NUMERICAL DATA FOR MIMIC PROGRAMS 

Numerical values of variable names a programmer may have 
used to designate constant coefficients, initial conditions, 
arbitrary functions, etc., are supplied to the processor on 
data cards following the END statement of the program. 

There are 4 types of data used in MIMIC programs and 
information concerning each type must appear in 2 places in 
a job. First, the name of the data must appear in a CON, 
PAR, CFN, or PFN statement and the corresponding value must 
appear on a data card. 

Values for data declared by a CON (Constant) statement 
stay fixed for all runs. One data card must appear in the 
data section for each CON statement in the program section. 

Values for data declared by a PAR (parameter) statement 
may change from run to run. One data card must appear for 
each PAR statement in the program section for each run de- 
sired. 

Up to 6 variable names may be used as arguments of each 
PAR or CON statement. 

Data for arbitrary functions of 1 (or 2) independent 
variable (s) declared by CFN (constant function) or PFN 
(parameter function) statements are punched as 1 pair (or 
triple) of points per data card in the following format: 

- functions of 1 variable; independent variable punched 
in columns 1-12, dependent variable punched in columns 
13-24 and the data cards ordered in increasing value 
of the independent variable. 
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- functions of 2 variables; 1st independent variable 
punched in columns 1-12, 2nd indpendent variable 
punched in columns 13-24, dependent variable punch- 
ed in columns 25-36. The cards are arranged so that 
the values of the 1st independent variable form a non- 
decreasing sequence and at each value of the 1st in- 
dependent variable , the values of the 2nd independent 
variable must form an algebraically increasing se- 
quence. 

In addition the array must consist of at least 
2 values of the second variable for each value of 
the first and the array must contain points less 
than and greater than the expected values of both 
independent variables. 

The number of data cards for both is specified as the 
argument of the CFN or FFN statement. 

Finally, any CON and CFN statements must come before 
any PAR and PFN statements and these cards must be the 
first cards in a program. The order of the data cards 
corresponding to each statement of each type is the same 
as the appearance of the statements in the program. 



For example, if 



F 
G 



CON(A,B) 
CFN (3.) 
PFN ( 4 . ) 
PAR(X) 



are statements in a program, data cards might be 



12 



24 



36 



-4. 


10. 




card 1: 


constant values of A,B 


-1.0 


-2. 




2j 


specify constant 


0. 


0. 




3[ 


function F of 1 


1. 


2. 




4) 


variable 


$MIMC,DCH3 










2 










2. 


-2. 


17. 


H 




2. 


5. 


10.8 


e\ 


specify function G of 



14. 



3. 


-2. 


25.0 


4. 


6.5 


100.63 


-25.36 






1.5 


-3.0 


25.0 


2.5 


-2.0 


14.0 


2.5 


4.0 


73.0 


3. 


0.0 


28.5 


-50.2 







7 
8 
9 

111 
12 
13 
14 



2 variables for 1st run 

value of X for 1st run 

function G for 
2nd run 

value of X for 2nd run 



NOTE: If *MIMIC load option is used, it must follow constant 
function followed by number of parameter sets, in this 
case 2, in 14 format columns 1-4. 

VII. OUTPUT 

1. Description of Printed Output 

of an OUT statement every DT units of T from T = until a 
FIN statement terminates the run. More than one OUT statement 
is allowed per program. 

e.g. 0UT(A,X) will print columns of values of A and X. 

A value for DT may be 

(i) read in by a CON or PAR card 

(ii) calculated by operations in a MIMIC program 

(iii) assumed to be .1 if not specified by (i) or (ii) . 

Blank columns can be introduced by using a comma with 
no preceding argument. For example OUT(A,,X). 

Blank lines may be inserted in the output by an OUT 
statement with no arguments. That is OUT ONLY . 

Up to 6 column headings to appear at the top of the out- 
put may be specified by HDR statements. The actual characters 
to be printed are the arguments of HDR and can consist of 1-6 
alphabetic or numeric characters. Blank column headings can 



15 



be obtained by using a comma with no preceding argument and 
blank lines in headings can be obtained by using HDR alone. 

In addition the processor prints out the names and values 
oj. any data it reads in for eacn run . 

1. Description of Plotted Output 

The plotter will supply a plotted output with the 
independent variable as the abscissa and the function as the 
ordinate. 

Up to four functions can be plotted simultaneously from 
the data points which are stored in the plotting subroutine. 
300 data points can be stored for each function; however, if 
this number is exceeded a diagnostic will result. 

The Y axis has a length of 20.0 inches, therefore, the 
user must scale his output to fit this. A mean value of 
10.0 is desirable. 

The X axis is 30,0 inches long and it is possible to 
scale the X-array by a multiplier supplied by the user. The 
time scale is normally marked at 1 second per inch. 

The user has the choice of plotting :- 

1. Data points only— — PLO(0. ,SCL,--) 

2. Straight lines between points PLO(l.,SCL, — ) 

3. A quadratic curve fit between points — PLO(2.,SCL, — ) 

Care should be taken using the third option if dis- 
continuities are expected. 

The second parameter in the call to plot, SCL, is the 
multiplier used to expand the graph in the X direction and 
is normally 1. The annotations of the X axis are modified 
accordingly. It is possible to reduce the scale of the whole 
graph. This is sometimes useful for photographic purposes. 
The decimal fraction of the first parameter is used to divide 
all dimensions in the plot. 
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Thus PLO(2.2,- ) will give a quadratic curve fit with 

the graph reduced by 2 . If the decimal fraction is blank or 
zero, the graph is of normal size (20" x 30") . 

VIII. SOME SPECIAL MIMIC FEATURES 

1. Logical Control Variables 

In addition to variables which can take on numeric 
values in the limits |l.7| x 10 +38, MIMIC allows logical 
variables which can take on only the values TRUE or FALSE. 
These variables can be used to provide control over the eval- 
uation of expressions. 

Logical variables may be assigned values by the use of 
element FWS (function switch) and the value can be changed 
by any of the logic functions (see parts 3, 7 of Table 1). 

If the name of a logical variable is entered in the LCV 
field of a MIMIC statement, the expression on that card will 
be evaluated only when that logical variable has the value 
TRUE. That is, logical control variables may be used to 
select, under programmed control, which of several paths 
should be evaluated. 

For example, define a logical variable which is TRUE 
only at T = O and use it to control precomputation of con- 
stants, initial conditions, etc. (see example 2 , Appendix I) 

NOTE : If a logical control variable is used on statements 
containing expressions involving INT or FTR, these 
functions will not be initialized unless the variable 
is TRUE at T = 0. 

2. Integrator Mode Control 

The mode of each integrator in MIMIC can be indi- 
vidually set to one of reset, operate or hold. In reset 
mode, the output of the integrator is set equal to the ini- 
tial condition argument and no integration takes place. In 
operate, the integrator is updated by the integration routine. 
In hold, the input to the integrator is removed. 
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Two extra arguments which must have logical values are 
used with INT to control the mode. If these are not present, 
the mode is assumed to be "operate". 



For example 



X INT(Y,XO,C,D) 



is a completely specified integrator. Logical variables C,D 
control the mode of the integrator by the following scheme:- 



D \ 


TRUE 


FALSE 


TRUE 


OPERATE 


HOLD 


FALSE 


RESET 


OPERATE 



3. Integrator Limiting 

MIMIC contains a function LIM to limit between some 
lower and upper bounds the values that a variable may assume 
However, if the variable to be limited is the result of an 
integration, it is necessary to set the derivative of the 
variable to zero whenever the variable is equal to a limit 
value by use of function LIN before using it as input to an 
integrator. 

For example, if x is the output of an integrator but is 
to be limited such that L <_ x <_ U, MIMIC statements might be 

XDOTLM LIN(XDOT,X ,L,U) 

X INT (XDOTLM, XO) 



or simply 



X 



INT(LIN(XDOT,X,L,U) ,X0) 



A block diagram for this is 
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x 



-» 



X vp 






LIN 



x (limited) 



x, 






X 



4. Implicit Functions 

MIMIC contains an element IMP for solving equations 
of the form x = f (x) . An iterative method is used and thus 
a starting value must be given by a PAR statement (or assigned 
in a calculation performed before any integration starts) . 



For example 



X 



PAR(X) 
IMP(X,SIN(X)+Y) 



would solve x = sin x + y where y is specified elsewhere in 
the program. 

Some computational efficiency can be obtained with im- 
plicit functions by using dummy variables since the processor 
function-sorting routine will insure that all computations 
involving the result of an IMP will be placed after the IMP. 



eg. 



PAR (XI) 



X IMP (XI, SIN (XI) +Y) 



NOTE: x = f (x) is solved by the scheme 
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x n = supplied value 



x 



n+1 



c = 

n 



f (x ) -ex 
n n n 



1 - c 



n 



f (x ) - f (x .) 
n n-1 



x - x n 
n n-1 



until 



x - f (x ) 
n n 



x 



n 



<.5 x 10 



-5 



•d . Time Delay 

The expression TDL(A,B,C) has the effect of delaying 
the variable A by B units of T. Here, B may be any expression; 
C, if specified must be a literal or constant, and represents 
the number of values of the input variable A stored at any 
one time for providing the delay. If C is not specified, it 
is assumed to be 100. 

If T <_ B, the result of TDL is the value of A at T = 0. 

6 . MIMIC Subprograms 

MIMIC language subprograms present a way of extend- 
ing the list of available functions for a given program. In 
particular, it allows functions with multiple outputs to be 
generated. 

A subprogram is defined as those statements occurring 
between BSP (begin subprogram) card and an ESP (end sub- 
program) card. The subprogram name is entered in the result 
field of the BSP and ESP cards. The dummy inputs to the sub- 
program are specified as arguments of the BSP element and the 
dummy outputs as arguments of the ESP element. Up to 6 
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inputs and outputs may be used and must be single variable 
names . 

A subprogram is used by entering its name in the result 
field of a CSP (call subprogram) card with the actual inputs 
listed as arguments of the CSP element and the actual outputs 
as arguments of the RSP element. Actual outputs must be 
single variable names but actual inputs may be expressions. 

The following rules apply in writing subprograms: 

(i) A dummy variable must be defined as an input or 
a result before it can be used as an argument. 

(ii) The names used in a subprogram must not be used 
elsewhere in the program. 

(iii) The RSP card must be the card immediately after 
the CSP card. 

For example, we define a subprogram SAMPLE to evaluate 

x, = x cos z + y sin z 

and y, = -x sin z + y cos z. 

The inputs are x,y,z and the outputs are x, ,y n . 

SAMPLE BSP(X,Y,Z) 
XI X*COS(Z) + Y*SIN(Z) 
XI -X*SIN(Z) + Y*COS(Z) 
SAMPLE ESP(X1,Y1) 

We could use SAMPLE in several places in a program pos- 
sibly as follows: 

SAMPLE CSP(U1,U2,U3) 
RSP(V1,V2) 



SAMPLE CSP(2.*T,45. ,PI) 
RSP(P1,P2) 
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Restrictions ; 

(i) The functions FTR and INT may not be used in 
subprograms . 

(ii) Arbitrary functions defined by CFN or PFN 
statements cannot be used as input to a 
subprogram; however, they may be used in 
FUN expressions within the subprogram. 

IX. CARD FORMATS 

1 . Program Cards 



I 2 



10 



15 



19 



72 



LCV 
field 



Result 
field 



Expression field 



I 



Any symbol in column 1 results in the entire card being 
treated as a comment. 

A single variable name only can appear in the LCV or 
Result field. 

2 . Data Cards 



12 



24 



36 



48 



60 



1st data 
field 



2nd data 
field 



72 



6th data 
field 



Up to six 12-column fields allowed per card. 



For CON and PAR statements, the value corresponding to 
the jth argument must be punched in the jth data field of the 
data card (1 <_ j <_ 6) . 

Data is punched in FORTRAN floating point format E12.0, i.e. 

(i) possible prefixed minus sign. 

(ii) explicit decimal point. 
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(iii) possible postfixed exponent of forms E+dd, E+d, 
Edd, Ed where d is a digit. 

(iv) no more than 12 characters per number. 

7 
eg. -1.72E7 means -1.72 x 10 

Note : An E-style number must be right justified in its data 
field. 

For example, for CON(A,B,C,D) a data card might be 



12 



24 



36 



48 



24.0 



-.0005 



■1.25E+10 



0.0 



Only the following errors will be detected and identified 
by the processor and will cause deletion of the execution 
phase. 

(i) Variable names which have not been assigned a value, 

(ii) Algebraic loops, 

(iii) A BSP card and no ESP card, 

(iv) A CSP card not followed immediately by an RSP card. 

(v) Use of a function not contained in Table 1. 

In addition, the usual execution time F-IV error mess- 
ages might appear. 

XI. SPECIAL POINTS AND SUGGESTIONS 

1. Integration 

(i) A maximum of 95 integrators is allowed. 

(ii) If DTMAX, DTMIN, DT are not specified in the 
program, the following values are assumed: 
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DT = DTMAX =0.1 

DTMIN =0. 
(iii) The processor uses a variable step, high order 
Runge-Kutta routine for integration. The step 
size may range from DTMIN to DTMAX. Setting 
DTMIN equal to DTMAX results in fixed step 
length integration. 

(iv) Use of hybrid functions FLF, MMV, TAS, ZOH re- 
quire that DTMAX = DTMIN. 

(v) The integration routine will automatically vary 
the step size in order to satisfy an error cri- 
terion. Thus, problems involving discontinuities 
generated by function switches, limiters, ran- 
dom number generators, etc. can cause trouble 
as the routine could continue to reduce the 
step size indefinitely, if a positive DTMIN is 
not specified, in order to integrate across 
discontinuities . 

It is impossible to give a rule for select- 
ing a DTMIN but the following approach is sug- 
gested. 

-4 
Try DTMIN = 10 DT and compare the results 

with another run with DTMIN one tenth the first 
value. If the results agree well enough, neither 
value inhibited the integration routine to any 
appreciable extent. If the results differ sig- 
nificantly, decrease DTMIN and rerun. Always 
try to use the largest value of DTMIN consistent 
with satisfactory results in order to decrease 
computer time . 

2 . Machine Fault Conditions 

A diagnostic of the following form will be returned every 
time a machine fault occurs as follows:- 
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ERROR IN MIMIC CALLED FROM XXXXX YYYYY 
Where XXXXX is a machine address 

and YYYYY is a. exponent overflow, 

b. divide fault, or 

c. arithmetic overflow. 

3. Control of Printed Output 

The processor outputs 60 lines per page on the printer. 

Thus if a DT = 10 were selected for integration over a 

5 
range as small as < T < 6 about 10 pages of output would 

be printed. It might be possible that only the final values 

for T = 6 are of interest. The following are suggestions 

for controlling the volume of output: 

(i) Specify DTMAX to control integration and specify 

(ii) Use logical control variables in the LCV field 
of OUT statements to produce selected output. 

For example: 



CONOUT FSW(T-B, FALSE, TRUE, TRUE) 



CONOUT OUT(T,X) 



will produce output only when T > B. 

(iii) Use more than one FIN statement to terminate 
the run. For example, terminate if some con- 
dition, say stable values of a variable, has 
been achieved or total time has elapsed. 
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4. Errors in Programs 

Section X details the only source language errors de- 
tected by the processor. Thus a programmer should take care 
in keypunching since the processor will evaluate expressions 
which are syntactically incorrect or outright meaningless. 
For example the processor will produce a value for 2.*TRUE. 

In compiling a program, the processor reduces all com- 
plex expressions to a sequence of calculations using the basic 
elements of Table I and lists this simplified program below 
the input program. A perusal of this simplified listing will 
often detect errors in the original program. 

XII. SAMPLE CONTROL CARDS FROM BATCH 

7 

8 JOB, < JOB #>,<USER £>,<LISTING ID> 

7 

8 TIME = <time limit> 

g MFBLKS = <File Block limit> 

7 

g *MIMIC,< PARAMETER STRING> 

MIMIC Program 
Constants 
Constant Functions 



Data card - number of parameter sets (F0RMAT(I4)) 

Parameter functions 

Parameters 



1 LOGOFF 



Note if plotted output is required, EQUIP a unit 
to a plotter; label it and have a PLOT = LUN PARA- 
METER on the control card. 
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g JOB, < JOB #>, <USER #> , <LISTING ID> 



7 



TIME = <TIME LIMIT> 
MFBLKS = <FILE BLOCK LIMIT > 
*MIMIC, <PARAMETER STRING> 

MIMIC PROGRAM 
CONSTANTS 
CONSTANT FUNCTIONS 

*MIMICLOAD, <PARAMETER STRING> 

DATA CARD - Number of parameter sets (FORMAT 14) 
PARAMETERS 



g LOGOFF 
XIII. SAMPLE CONTROL STATEMENTS FROM TELETYPE 



#<JOB#>,<USER#> 

#TIME=<TIME LIMIT> 

#MFBLKS=<FILE BLOCK limit> 

#* MI MIC, I =NAME,< PARAMETER STRING> 

#LOGOFF 

XIV. CONTROL CARD PARAMETERS 

Options on the call to the compiler are separated by 
commas and each may be any number of letters long, and may 
be followed by an =XX where XX is some logical unit number. 
For *MIMIC the options currently available are: 

ASSEMBLY produces an assembly listing on the requested 
unit. If no unit is specified, the assembly 
listing appears on the same unit as the source 
listing (Logical unit 61 if neither the LIST 
or ASSEMBLY options have specified a unit) . 
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CONSTANT 



FUNCTION 



INPUT 



LIST 



PLOT 



OUTPUT 



produces a listing of the constants defined on 
CON or CFN cards. If no unit is specified, 
listing defaults to other units in the same 
order as the FUNCTION listing. 

produces the SORTED FUNCTION LISTING GENERATED 
by the compilier. If no unit is specified use 
the assembly unit if specified or default to 
the list unit if specified or USE 61. 

specifies the unit or file from which the 
source program is to be read. (Logical unit 
60 is assumed if no unit is specified or if 
the option is not present.) The compiler will 
automatically recognize and unblock EDIT and 
COSY compressed records in addition to stan- 
dard Hollerith card images. As an additional 
feature, the name of a saved file may be used 
in place of a logical unit number on this op- 
tion. MIMIC will rewind the input unit before 
reading if a name was specified as the input. 

specifies that a listing of the source is to 
be written out on the specified unit (Logical 
unit 61 if no unit was specified.) 

produces the plot as specified on PLO cards. 
(If PLO cards are present and the plot option 
is not used the plot will be suppressed.) 
Assumes Logical unit 55 if none specified. 
(However it is a good idea to equip a unit to 
a plotter before starting a run and labeling 
the unit.) 

the unit to list the output generated by the 
source program. Logical unit 61 is assumed 
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run indicates a unit (Logical unit 56 if none is 

specified) on which binary output is to be 
generated. The RUN unit is released initially 
by the compiler, and on compiling the last 
program element, is automatically rewound and 
loaded and execution is initiated if no seri- 
ous errors were encountered in compilation. 

XCT indicates a unit (Logical unit 56 if none is 

specified) on which binary output is to be 
generated. MIMIC releases this unit before 
compiling. 

If the RUN option is not present and either the plot or 
output option is present they are ignored. If the source 
program input and the data for the source program are on dif- 
ferent Logical units the RUN cannot be used. 
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For *MIMICLOAD the options currently available are: 

INPUT specifies the unit or file from which the data 

for the source program is to be read. (Logical 
unit 60 is assumed if no unit is specified or 
if the option is not present.) The compiler 
will automatically recognize and unblock EDIT 
and COSY compressed records in addition to 
standard Hollerith card images. As an addi- 
tional feature, the name of a saved file may 
be used in place of a logical unit number on 
this option. The compiler will rewind the in- 
put unit before reading. 

OUTPUT the unit to list the output generated by the 
source program. Logical unit 61 is assumed 
if not specified. 

PLOT produces the plot as specified on PLO cards. 

(If PLO cards are present and the plot option 
is not used the plot will be suppressed.) 
Assumes Logical unit 55 if none specified. 
(However it is a good idea to equip a unit 
to a plotter before starting a run and label- 
ing the unit.) 

XCT unit on which to find binary output from 

*MIMIC. If no unit specified assumes Logical 
unit 56. 

NOTE: *MIMICLOAD is a loader that loads binary output pro- 
duced by *MIMIC and then starts execution. 

SAMPLE COMPILER CALLS 

*MIMIC,R,L 
*MIMIC,X,I=10,L,F=5 
*MIMICLOAD,X, J=l , 0=45 ,P=13 
*MIMIC,R,L,A,F,C,P 
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APPENDIX I 



Example 1 : The equations for a bouncing ball dropped from 
rest at height h are x = -g, x (0) = 0, x(0) = h where x (t) 
represents the height at time t. But whenever x - 0, the 
velocity is changed in direction and decreases in magnitude 
by proportion e, the coefficient of elasticity. 

A possible job is 



1 2 



10 



19 



24 







| 


CON (TS 
PAR (E, 


rOP, DTMIN) 
H) 








ZEROHT 


T-i ri t.t / \r rn 

row \a,x 






ZEROHT BECOMES 


TRUE WHEN 


X REACHES 


ZERO AND IS FALSE FOR X POSITIVE 


! 


XD0T 


INT(-32 


.2,0.) 




ZEROHT 


XDOT 


-E*XDOT 










X 




INT(XDO 
HDR (TI 
HDR 

OUT (T, 
FIN (T, 
END 


r,H) 

yiE, SPEED, HEIGHT) 

XDOT,X) 
rSTOP) 




10.0 
j 






1.0E-4 






0.09 






.15E+02 








.8 




>0.0 








1 


.5 




>5.0 









The use of logical variable ZEROHT ensures that we 
account for elasticity and direction only at times when 
x < 0. ( < used should the integration routine slightly 
overshoot to negative x) 
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Example 2 : Illustration of subprograms and precomputation 
of constants. 

Evaluate the complex-valued variable 

z = P(s)/q 

2 
where P(s) = us - vs + w = (us-v)s + w 

ia 
u = e 

v = log(sinB) 

w = acos3 + iftsina 

q = a + i3t 

s = t(a + i$) 

for t=0 (.2) 2. (a, 3 are real numbers) 

Since u, v, w remain constant for all t, we will compute 
them onlv once when t = bv usincr a loaical control variable 

~ -t r -*■ -* •* 

which has the value TRUE only when T = 0. 

Subroutines COMPAS, COMPMD will be defined as follows: 
COMPAS - number of inputs = 5 

- inputs 1, 2, 3, 4 are real and imaginary parts of 
2 complex numbers 

- input 5 is a logical-valued argument 

- if its value is TRUE, COMPAS adds 

- if its value is FALSE, COMPAS subtracts 

- number of outputs = 2 

- the real and imaginary parts of the result. 
COMPMD - number of inputs = 5 

- inputs 1, 2, 3, 4 are real and imaginary parts 
of 2 complex numbers 
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input 5 is a logical-valued argument 

- if its value is TRUE, COMPMD multiplies 

- if its value is FALSE, COMPMD divides 



- number of outputs = 3 



The program: 



- outputs 1, 2 are real and imaginary parts of 
the result 

- output 3 - a logical-valued variable 

- its value is FALSE only if the input divisor 
was zero for a complex division, (Outputs 
1, 2 = in this case) Otherwise it has 
the value TRUE. 



1 2 



10 



19 



EEFINE 



DEFINE 



SUBPROGRAM 

COMPAS 

SIGN 

RC 

IC 

COMPAS 
SUBPROGRAM 

COMPMD 

SIGN1 

R3 

13 

NOTL 

LOUT 

ZERO 



CON (ALPHA, BETA) 

COMPAS 
BSP (RA,IA,RB,IB,LIN) 
LSW(LIN,1. ,-1.) 
RA+SIGN*RB 
IA+SIGN*IB 
ESP (RC,IC) 

COMPMD 
BSP (R1,I1,R2,I2,L0GIN) 
LSW (LOGIN, 1. ,-1.) 
Rl*R2-SIGN*Il*I2 
SIGN1*R1*I2+I1*R2 
NOT (LOGIN) 
TRUE 
FALSE 
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1 2 



10 



19 



r~nr 



NOTL 
NOTL 

ZERO 

i 

JZERO 
JNOTL 



^RECOMPUTE 

I 

i 

ZEROT 
JZEROT 
jZEROT 
1ZEROT 
IZEROT 
JZEROT 
JZEROT 



"denom 

ZERO 

*3 

13 

LOUT 

COMPMD 

U,V,W, 
ZEROT 
RU 
IU 
RV 
IV 
RW 
IW 
PT 



CALCULATION OF Z 



IOGVAL 



COMPMD 



COMPAS 



COMPMD 



JCOMPAS 



COMPMD 



\ 



R2*R2+I2*I2 

FSW ( DENOM , TRUE , FALSE , TRUE ) 

R3/DENOM 

I 3 /DENOM 

ZERO 

ESP (R3, 13, LOUT) 

; DT 

FSW (T , TRUE , TRUE , FALSE ) 

JCOS (ALPHA) 

SIN (ALPHA) 

LOG (SIN (BETA) ) 

o. 

ALPHA*COS (BETA) 
BETA*SIN (ALPHA) 
0.2 

CSP (RU,IU,T*ALPHA,T*BETA, TRUE) 

RSP(RP,IP,LOGVAL) 

CSP (RP, IP, RV, IV, FALSE) 

RSP(RQ,IQ) 

CSP ( RQ , IQ,T* ALPHA ,T*BETA, TRUE) 

RSP(RS,IS,LOGVAL) 

CSP(RS,IS,RW,IW,TRUE) 

RSP(RT,IT) 

CSP (RT , IT , ALPHA*T , BETA*T , FALSE) 

RSP(RZ,IZ,LOGVAL) 

OUT(T,RZ,IZ) 
FIN(T,2.) 

Iend 



LOGVAL is used as an LCV on the OUT statement so no results 
are printed in case a zero division was attempted in COMPMD. 
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APPENDIX II 

Solution of a Partial Differential Equation by Finite Approximations 

The equation for the concentration C (r,t) of a liquid 
component diffusing radially out of a solid particle cylinder 
of radius R is 



!£= D 

9t 



8 2 C 1 8C 

2 

3r r 8r 



D is the diffusivity. 

There is an initial uniform concentration distribution 
C* throughout the cylinder and at time t = we allow the 
material to commence diffusing by putting C(R,0) to zero and 
maintaining it at zero. 

The analytic solution to the problem is 

C(r,t) = ? 2C * J n (a. r ) exp (- i t) II-l 

u x 2 

i=l a. J, (a. ) R R 

(See Reddick and Miller: Advanced Mathematics for Engineers, 
J. Wiley, p. 319, Eq. 21) 

a. are the positive zeros of J Q (x) 

To acquire numerical values of C(r,t) using MIMIC we make 

the approximations 

9 2 
C C-C °CC2C+C 
n n+1 n-1 , n = n+1 n n-1 

2 2 

9r 2Ar dr (Ar) 

T~. Contributed by Professor T. Fahidy, University of Waterloo 
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where C = C(nAr,t). The partial differential equation is 
turned into a set containing n ordinary differential equations . 
This is a conventional technique for analogue computers. 



Now we have 
dC 



n D 



dt (Ar) 2 



(1 + — ) C , n - 2C + (1 — ) C . 

n n+1 n n n-1 
2n 2n 



XX "~~ If Z^ Jjt 7^ •«• 

for the concentration on concentric circles of radius nAr. 

Taking C* = 1, R = 0.5, r= 0.1, n = 5, D = 10~ 5 we get 
C. = 10" 3 (1.5C - 2C + 0.5 C„) 

J. Z J. u 

C 2 = 10~ 3 (1.25C 3 - 2C 2 + 0.75 C^) 

C 3 = 10" 3 (1.167C 4 - 2C 3 .+ 0.833 C 2 > 

C 4 = 10~ 3 (1.125C 5 -- 2C 4 + 0.875 C 3 ) (C 5 = = boundary value) 

A MIMIC program to integrate this system of equations for 
C. in t = (10.) 200. along with all computer output follows. 

Total 3300 time for this job was 1 minute 21 seconds. 
The value of C(0.4, 200) calculated from II-l is 
0.872488 in contrast with the MIMIC - computed value of 
0.813412. The error (about 7% may be further reduced 
by considering smaller increments of Ar. This is a case 
where accuracy and program complexity have to be bal- 
anced. 
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MIMIC SOURCE -LANGUAGE PROGRAM 



The Radial Diffusion Problem with 5 Increments 

C0N(1A,2A,3A,2N) 

C0N(1B,2B,3B,4B) 

CON(DT,DTMAX) if 
1C INT(1A*2C-2N*1C+1B',1. ) 
2C INT(2A*3C-2N*2C+2B*1C,1.) 
3C INT(3A*4C-2N*3C+3B*2C,1.) 
4C INT(-2N*4C+4B*3C,1.) 

FIN(T,200.) 

HDR(T,C1,C2,C3) 

HDR( f C4) 

HDR 

0UT(T,1C,2C,3C) 

0UT(,4C) 

END 
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***FUNCT I ON -LANGUAGE PROGRAM GENERATED*** 



1J!'H i-i<~V 

1 

2 

3 

4 

5 

6 

7 

8 

9 
10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 



T-\T-| <"* TTT IT! 


i-nrpvT 
r ilM 


TV 


B 






CON 


1A 


2A 


3A 2N 




CON 


IB 


2B 


3B 4B 




CON 


DT 


DTMAX 




103 


MPY 


1A 


2C 




104 


MPY 


2N 


1C 




105 


SUB 


103 


104 




106 


ADD 


105 


IB 




1C 


INT 


106 


1. 




108 


MPY 


2A 


3C 




109 


MPY 


2N 


2C 




110 


MPY 


2B 


1C 




111 


SUB 


108 


109 




112 


ADD 


111 


110 




2C 


INT 


112 


1. 




114 


MPY 


3A 


4C 




115 


MPY 


2N 


3C 




116 


MPY 


3B 


2C 




117 


SUB 


114 


115 




118 


ADD 


117 


116 




3'C 


INT 


118 


1. 




120 


MPY 


2N 


4C 




121 


MPY 


4B 


3C 




00 


NEG 


120 






122 


ADD 


00 


121 




4C 


INT 


122 


1. 






FIN 


T 


200. 






HDR 


T 


CI 


C2 C3 




HDR 




C4 






HDR 










OUT 


T 


1C 


2C 3C 




OUT 




4C 






END 









***FURTHER DIAGNOSTICS AND EXECUTION FOLLOW*** 

1A 2A 3A 2N 

1.50000E-03 1.25000E-03 1.16700E-03 2.00000E-03 

IB 2B 3B 4B 

5.00000E-04 7.50000E-04 8.33000E-04 8.75000E-04 

DT DTMAX 

1.00000E 01 1.G0O00E-01 
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T 




CI 
C4 




C2 






C3 




U . 






1. OOOOOE 00 
1. OOOOOE 00 




.OOOOOE 




1. 


.OOOOOE 


00 


1. 


.OOOOOE 


01 


9.99998E-01 
9.88860E-01 


9. 


.99998E- 


01 


9. 


.99934E- 


-01 


2. 


.OOOOOE 


01 


9.99997E-01 
9.77940E-01 


9. 


.99995E- 


-01 


9. 


.99742E- 


-01 


3. 


.OOOOOE 


01 


9.99994E-01 
9.67234E-01 


9. 


.99989E- 


-01 


9. 


.99429E- 


-01 


4. 


.OOOOOE 


01 


9.99992E-01 
9.56737E-01 


9, 


.99978E- 


-01 


9. 


.98999E- 


-01 


5. 


.OOOOOE 


01 


9.99990E-01 
9.46443E-01 


9, 


.99962E- 


-01 


9, 


.98458E- 


-01 


6. 


.OOOOOE 


01 


9.99987E-01 
9.36348E-01 


9, 


.99938E- 


-01 


9. 


.97810E- 


-01 


"7 


.OOOOOE 


01 


9.99986E-01 
9.26447E-01 


9, 


.99907E= 


=01 


9, 


,97058E= 


-01 


8. 


.OOOOOE 


01 


9.99984E-01 
9.16735E-01 


9, 


.99866E- 


-01 


9, 


.96208E- 


-01 


9. 


.OOOOOE 


01 


9.99980E-01 
9.07208E-01 


9. 


.99814E- 


-01 


9. 


.95263E- 


-01 


1. 


.OOOOOE 


02 


9.99976E-01 
8.97860E-01 


9. 


.99752E- 


-01 


9, 


.94228E- 


-01 


1. 


.lOOOOE 


02 


9.99972E-01 
8.88689E-01 


9, 


.99677E- 


-01 


9. 


.93105E- 


-01 


1. 


.20000E 


02 


9.99965E-01 
8.79688E-01 


9, 


.99590E- 


-01 


9, 


.91899E- 


-01 


1. 


.30000E 


02 


9.99958E-01 
8.70856E-01 


9. 


.99488E- 


-01 


9. 


.90613E- 


-01 


1. 


.40000E 


02 


9.99949E-01 
8.62186E-01 


9, 


.99372E- 


-01 


9. 


.89250E- 


-01 


1, 


.50000E 


02 


9.99939E-01 
8.53677E-01 


9. 


.99241E- 


-01 


9, 


.87815E- 


-01 


1. 


.60000E 


02 


9.99927E-01 
8.45323E-01 


9. 


.99095E- 


-01 


9. 


.86309E- 


-01 


1. 


.70000E 


02 


9.99912E-01 
8.37121E-01 


9. 


.98932E- 


-01 


9, 


.84736E- 


-01 


1. 


.80000E 


02 


9.99896E-01 
8.29067E-01 


9. 


.98752E- 


-01 


9. 


.83099E- 


-01 


1. 


.90000E 


02 


9.99877E-01 
8.21159E-01 


9. 


.98555E- 


-01 


9. 


.81400E- 


-01 


2. 


.OOOOOE 


02 


9.99855E-01 
8.13392E-01 


9, 


.98341E- 


-01 


9. 


.79643E- 


-01 
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APPENDIX III 



Example of Integrator Mode Control and Hybrid Functions 

This problem is designed to illustrate the use of in- 
dividual mode control of the integrators and the use of cer- 
tain other logical elements in the solution of problems 
which more conventionally might be solved on a hybrid analogue- 
digital computer installation. The problem chosen is that 
of hitting a target at a known distance with a projectile 
fired at a certain muzzle velocity and at some angle of ele- 
vation to the horizon. 




TARGET X */ 
ERROR" 



Actual dynamic equations describing such a system are 
quite complex when all the aerodynamic forces on the pro- 
jectile are considered. The method of solution illustrated 
by this problem is general enough to be used for the deter- 
mination of the firing parameters in complex non-linear rep- 
resentations of the projectile equations. However, for the 
sake of computer time and so as not to obscure the method 
of solution by using an unnecessarily difficult equation set, 
the actual equations chosen have been simplified to the point 
of triviality. 



1. Contributed by E. M. Hinchley 
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The equations (neglecting all forces except gravity) are: 



x = V 

y 



cose 





-g 



x(0) 


= 


y(0) 


= V 


y(0) 


= 



sin( 



The desired point of impact is designated by x = RANGE 
and y is the variable representing the height of the pro- 
jectile. V is the muzzle velocity, and 6 is the angle of 


elevation (see sketch) . For simplicity V will be assumed 

constant and 9 will be the parameter we must determine. 
The equations of the projectile are so simplified that they 
are easily solved by direct integration. However, we are 
attempting to illustrate a more general method of solution- 
an iterative analog computer-like solution. 

The equations are programmed in MIMIC as they would be 
for an analog computer solution. A value of 8 is assumed and 
the analog solution is operated until y <_0(and x>0) . The 
value of x at y=0 represents the distance of the point of 
impact. The analog solution is placed in the HOLD mode 
while the error function x-RANGE is calculated, and a new 
value of 9 equal to the previous 6 plus K (x-RANGE) is deter- 
mined. K is a carefully chosen constant which affects both 
the convergence properties and the stability of the solution. 
The analog solution is next RESET to the new initial condi- 
tions and then placed in OPERATE and the new point of impact 
is determined. The process is repeated until a point of 
sufficiently close to the desired RANGE is found. 

The modes of the integrators are controlled by logical varia- 
bles C and D according to this table: 



D Nv 


TRUE 


FALSE 


TRUE 


OPERATE 


HOLD 


FALSE 


RESET 


OPERATE 



Through the use of two function switches in the program, a 
logical variable Bl is obtained such that Bl is FALSE for y<^0, 
x > and TRUE otherwise. Another logical variable S becomes 
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FALSE at the same time. It was decided arbitrarily that the 
solution should HOLD for 1 second after impact, and then RESET 
for 1 second after that. The timing was obtained through the use 
of monostable multivibrators (MMV) . The logical variable H is 
TRUE when S is TRUE and becomes FALSE one second after S becomes 
FALSE. Logical variable becomes FALSE 2 seconds after S becomes 
FALSE. The integrator control variables C and D are obtained by 
"ANDing" and "ORing" combinations of the above variables to give 
the proper mode cycling from OPERATE to HOLD to RESET to OPERATE. 
The logical functions are: 

S = H.O-Bl+Bl.O" j • is AND 

C = S-O-H+S-H + is OR 

D = H^O I x is NOT x 

Initially Bl is TRUE, and S is made TRUE, is assuring that 
H and are TRUE so that C and D are TRUE and the solution will 
OPERATE. 

These logical functions were obtained through the use of a 
simple truth table. Some other simpler methods of controlling the 
integrator modes were tried, but since the 7040 is a serial mach- 
ine, it was found that the simulation of parallel logic required 
a great deal of care. 

The variable THETA ( ) is updated in the HOLD mode only 
through the use of a track and store device (TAS) and logical 
control variables. It was necessary to use equal integration 
intervals because of the logical devices being used. Despite the 
simplicity of the equations and the use of a fairly large integra- 
tion interval, solution times were surprisingly large on the 7040. 
Convergence towards a solution was obtained. A program listing 
including comment cards follows: 
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PARAMETERS 

EQUAL INTEGRATION 
DTMIN 
DTMAX 

LOGICAL VARIABLE 



PAR ( RANGE , VO , THETAO , G , K , DT ) 

INTERVAL 

EQL(.Ol) 



■pat f r\Ti\AT\l "\ 



AT THE END OF EACH ITERATION 



TQ TRUE 



LOGICAL 
TQ 
NTQ 



Bl BECOMES FALSE 

B FSW(Y, FALSE, FALSE, TRUE) 

Bl FSW(X, TRUE, TRUE, B) 

AT T = 0, CONTROLS PRECOMPUTATION OF CONSTANTS 

TQ FSW(T, FALSE, TRUE, FALSE) 

NTQ COM(TQ) 

VARIABLES FOR MODE CYCLING FOLLOW 

S EQL(TRUE) 

S IOR(AND(H,0,Bl) ,AND(B1, COM(O) ) ) 

H MMV(S,1. 0) 

MMV(S,2.0) 

C IOR(AND(S,0,H),AND(COM(S),COM(H))) 



NEXT 3 CARDS 
TH 
CP THETA 

TQ THETA 

ACTUAL PROJECTILE 

X 
INTEGRATE -G TWICE 



COM(C) 

ft TlTT-v /TT /-\ \ 



COMPUTE THETA, UPDATING IT IN THE HOLD MODE 
TAS ( THETA , C , THETAO ) 
EQL(TH+COR) 
EQL( THETAO) 

EQUATIONS ARE SOLVED BY THE NEXT 2 STATEMENTS 
INT(VO*COS(THETA) ,0.0,C,D) 
WITH RESPECT TO TIME 



USING MODE CONTROLS 

Y INT (INT (-G,VO*SIN (THETA), C,D) ,0. ,C ,D) 

ERROR FUNCTION AND A CORRECTION FOR THE NEXT VALUE OF THETA ARE FOUND 



FINISH 
CP 



ERR 
COR 
CONDITIONS 



SUB (X, RANGE) 
MPY(K,ERR) 



FIN(100.0,ABS(ERR)) 

FIN(T,200.) 
TITLES AND OUTPUT CARDS 

HDR(T,X,Y, THETA) 

HDR 

OUT(T,X,Y, THETA) 
NUMERICAL VALUES FOR PARAMETERS FOLLOW END CARD 

END 
2000.0 750.0 0.8 32.2 



0.00001 



0.2 
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APPENDIX IV 

The Simulation of a Servo Control System 

Figure 1 shows the control system to be simulated using 
three settings of the compensation network bandwidth denoted as 
XI, Xj and Xj Q . 

The transfer function 

(1 + ST X ) (1 + ST 3 ) 



is reduced to 



T,T 



S(l + ST 2 ) 



T l + T 3 ~ T 1 T 3 - T 2 



13 , 1 _,_ T 
+ S + 2_ 



T 2 1 +'T,S. 

and the equivalent mimic elemental representation of the servo 
system is shown in Figure 2. Use is made of the first order 
transfer function, FTR, the integrator INT and the usual arith- 
metic elements. Listing of the source program is shown in 
Figure 3 and the plotted output Figure 4. Total computer run 
time on the CDC 3300/MASTER was 1 minute 38 seconds. 







COMPENSATION 


CONTROLLER 

OUTPUT 

VOLTAGE 

/~V_ 


VELODYNE 
TRANSFER 
FUNCTION 




AXIS 

POSITION /^ 


T - 


K(l+ .T^Cl + sT 3 ) 


AXIS 
POSITION 


DEMAND +/ 


g J _ 




k 
S 






»( 1 + sT 2 ) 


^— 






»- 



PARAMETER 



EL 



45.4 



.222 



XI 



15 



165 



2.5 



AZ 



25 



.400 



X 1/3 



EL 



4.54 



.222 



AZ 



2.5 



.400 



47.4 



522 



7.91 



X 1/10 



EL 



.454 



,222 



AZ 



.25 



.400 



150 



1650 



25 



UNITS 



VOLTS/SECOND/DEGREE 



DEGREES/ SECOND /VOLT 



SECONDS 



SECONDS 



SECONDS 



FIGURE 1 - COMPENSATED TYPE 2 POSITION LOOP 






MIMIC EQUIVALENT OF FIGURE V 




T l + T 3- T lV T 2- T 2 



FIGURE 2 
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TO RUN FROM BATCH THE DECK LOOKS LIKE: 

000100 #<JOB#>,<USER#>,<USER ID > 

000200 ?EQUIP,l=PLOT 

000300 Label, 1/ <user id> 

000400 2#mimic,l,p=1,r 

000500 this is a simulation of a precision antenna 

000600 control system 

000700 CON(THETAI,DT) 

000800 PAR(T1,T2,T3,K1,K2) 

000900 Al SUB (THETAI, THETAO) 

001000 A2 MPY(A1,K1) 

001100 . Bl MPY.(A2.(T1*T3)/T2) 

001200 B2 INT(A2,0.) 

001300 B31 MPY(A2 5 T1+T3-T1*T3/T2-T2) 

001400 B3 FTR(B31,T2) 

001500 Dl ADD(B1,B2,B3) 

001600 D2 MPY(D1 5 K2) 

001700 THETAO INT(D2,0.) 

001800 FIN(T,120.) 

001900 HDR(TIME,RESP) 

002000 HDR 

002100 OUTCT, THETAO) 

002200 PLOC2.,. 25, THETAO) 

002 300 END 

002400 1.00000E 01 8.00000E-01 

002500 3 

002600 1.50000E 01 1.65000E 02 2.50000E 00 2.50000E 01 0.40000E 00 

002700 4.74000E 01 5.22000E 02 7.91000E 00 2.50000E 00.0.40000E 00 

002800 1.50000E 02 1.65000E 03 2.50000E 01 2.50000E-01 0.40000E 00 



FIGURE 3 



48 



" § 



ii 
<r 



it 

CO 



ii 
o 



II 
o 



I H- 



■ t ii ii i .i i «;^ '■> i i i i ■ '■ « • * 

xt-i oon «■«■« ° 30 '' 

S318diybA lN3QN3d3Q 




- g 



■ & 



ID 

Q_ 
I— 
ID 
O 



CE 
C£ 
CD 
O 

or 



